Spin-Strain Phase Diagram of Defective Graphene 
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Using calculations on defective graphene from first principles, we herein consider the dependence 
of the properties of the monovacancy of graphene under isotropic strain, with a particular focus 
on spin moments. At zero strain, the vacancy shows a spin moment of 1.5 \±b that increases to 
~2 fiB when the graphene is in tension. The changes are more dramatic under compression, in that 
the vacancy becomes non-magnetic when graphene is compressed more than 2%. This transition is 
linked to changes in the atomic structure that occurs around vacancies, and is associated with the 
formation of ripples. For compressions slightly greater than 3%, this rippling leads to the formation 
of a heavily reconstructed vacancy structure that consists of two deformed hexagons and pentagons. 
Our results suggest that any defect-induced magnetism that occurs in graphene can be controlled 
by applying a strain, or some other mechanical deformations. 



Since its discovery [lj, graphene has been shown to pos- 
sess a number of remarkable electronic and elastic prop- 
erties Q. Its electronic properties are mainly associated 
with the p z orbitals of carbon, while its elastic proper- 
ties are associated with the sp orbitals. Specifically, a 
Young's modulus of 1 TPa[3] has recently been measured 
for graphene, in agreement with the average Young' s 
modulus for carbon nanotubes of Y =1.25 TPa [J H|- 
Graphene can sustain elastic deformations as large as 
20%, and it is typically under a strain of several per- 
cent when deposited on surfaces [6] . This strain could be 
used to engineer the electronic properties of graphene. 
Unfortunately, a very high tensile force must be applied 
to defect-free graphene to modify the electronic struc- 
ture, because the spectrum of planar graphene remains 
metallic up to 10% stretching [7H9j. However, there are 
not many experimental conditions under which layers of 
graphene remain strictly planar. Due to the soft transver- 
sal phonons of 2D graphene, at certain temperatures the 
layer shrinks and ripples are formed. These ripples are 
present in free-standing layers of pristine graphene [lol [Tlj 
and for layers deposited on substrates. |12l, [L3| . The rip- 
pling of graphene induces midgap states [14[ in the elec- 
tronic structure that disappear when relaxations in the 
layers are carefully taken into account [l5j|. The relax- 
ations that accompany the rippling of the layer, and their 
influence on the electronic structure, are particularly im- 
portant in layers that are under a high degree of compres- 
sion. Strong rippling of this kind can also be induced by 
adsorbates [16| and defects jl7| such as that recently dis- 
covered for OH impurities [ll|, [l8|, [l9|. Understanding 
the interaction between defect-induced deformations and 
ripples is one of the main challenges now faced in the 
study of the electronic structure of graphene. 

In pristine graphene, the sp contributions are situ- 
ated at low energies, and it is the presence of defects 
that promotes them closer to the Fermi level. By the- 



oretical calculations, it has been shown that substitu- 
tional doping (20l - [23| and the presence of defects 0, HH 
in graphene, and in graphitic materials in general, pro- 
duce magnetism that is of interest in the potential use of 
these materials in spintronics. In experiments, the irradi- 
ation [26l-l29| and ion bombardment [30| of carbon-based 
materials create vacancies which indeed are linked with 
magnetic signals. The vacancies [24], [HI, [32| and edges 
[33L l34j present in graphene layers have been the focus 
of detailed theoretical studies. Although there have been 
both studies of the changes in the electronic structure 
induced by rippling and studies of defects in graphene, 
to our knowledge there have been no studies of the mag- 
netism of rippled graphene. Such a study is of particular 
interest when considered alongside the effect of strain on 
the structural and electronic properties of graphene. 

We herein combine the two perspectives and present 
a study of defective graphene under strain. We focus 
on a particular type of defect, namely related to carbon 
monovacancies, and find that these show a rich diversity 
of structural and spin phase-behavior as a function of 
strain. For zero strain, each vacancy has a solution with 
spin moment of 1.5 \±b- However, non-magnetic solutions 
become stable under a moderate compression of less than 
2% strain. This transition in magnetic moment may be 
traced to the modification in the local atomic structure of 
the defect as ripples develop in the layer. In calculations 
for compressions greater than 3%, the vacancy structure 
departs from its usual geometry under rippling, and this 
reconstruction ends with the formation of a different de- 
fect structure that has a central atom with strong sp 3 hy- 
bridization. All these structural changes are linked with 
the appearance of ripples as the layer is compressed, sim- 
ilar to those observed in recent experiments (35— 37] . We 
herein explore the relationship between the appearance of 
rippling patterns, the presence of vacancies, and changes 
in the electronic and magnetic properties of graphene. 
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We perform density functional calculations |4£ 
first principles using the SIESTA code 147 
use the generalized gradient app roximation [48j|, norm- 
conserving pseudopotentials [49| and a basis set of nu- 
merical atomic orbitals [47j. The integration over the 
Brillouin zone uses a well-converged Monkhorst-Pack k- 
sampling [50| that is equivalent to 70 x 70 points for the 
graphene unit cell. The structural optimizations use the 
conjugate gradient algorithm, which causes the forces to 
converge until they are lower than 0.05 eV/A. 




FIG. 1. (A) Graphene with vacancies under an isotropic com- 
pression of 1.2%. As an example we use the 10x10 unit cell 
(highlighted). The inset shows the geometry of the vacancies 
and the atomic labels. The local bending at the vacancies 
is described by the angle 612 between the pentagon and the 
plane defined by three C atoms around the vacancy, i.e. two 
equivalent atoms labeled 1 and a third atom labeled 2. Note 
the rippling of graphene sheets with vacancies at the saddle 
points. (B) A different vacancy structure for a compression 
slightly under 3%. It has two distorted hexagons and pen- 
tagons, while the central C atom shows strong sp 3 hybridiza- 
tion. 

Structure versus strain. While under tension the layer 
remains perfectly flat, and small compressions produce 
the spontaneous rippling of graphene with a character- 
istic deformation around the vacancy. [T] illustrates the 
rippled geometry after relaxing one of our models for a 
free-standing defective graphene layer under an isotropic 
compression of 1.2%. At zero strain, the monovacancy 
tends to undergo a J ahn- Teller-like distortion that low- 
ers its energy by ~200 meV: atoms of type 1 (see the inset 
of [Tj) reconstruct to form a pentagon with the neighbor- 
ing atoms, while atom 2 is left with the dangling bond 
responsible for the spin polarization. Although we herein 
consider a non-planar structure under strain, our findings 
at low compressions are similar to results previously ob- 
tained for flat graphene [38, 39|. However, atom 2 in our 



rippled structure is progressively lifted from the graphene 
surface by as much as ~1.0 A for strains slightly below 
3%, just before the occurrence of a strong vacancy recon- 
struction. The data of []J3 reveal a particularly striking 
result. In contrast to 2D vacancy, for which previous au- 
thors invariably described a 'pentagon-goggles' structure, 
our own calculations show that the vacancy is strongly 
reconstructed under rippling. For compressive strains 
greater than 3%, a different defect structure is produced 
altogether, consisting of two heavily distorted hexagons 
and pentagons and resembling the transition states of a 
planar vacancy movement. We shall now focus on the 
description of this region of transition for compression of 
up to 2.8%. 
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FIG. 2. (A) Maximum amplitude h(A) of rippling for 10x10 
graphene supercells versus strain in units of the equilib- 
rium lattice constant a . The filled circles denote defec- 
tive graphene, while the empty circles denote the pristine 
layer that has been recovered by adding a carbon atom to 
the vacancy. The insets show the corrugation patterns: the 
bright /yellow (dark/green) areas indicate the higher (lower) 
regions. (B) Bending angle 612 (see definition in[T]) as a func- 
tion of strain. The inset shows the behavior of the bending 
angles with size of supercell for different values of strain. The 
local induced deformation clearly depends on the presence of 
vacancies and on the size of the supercell. 

[2]A. shows the maximum amplitude of the rippling h as 
a function of strain for defective graphene (shown by the 
solid symbols). We removed the defects from each rippled 
structure by adding carbon atoms to the vacancy sites 
and relaxing the structure. We obtained rippled pristine 
graphene with maximum amplitudes as shown in [2] (open 
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symbols). The amplitudes and the topographic patterns 
(insets in [2JA.) show little difference between pristine or 
defective graphene, particularly for compressive strains 
greater than It seems that even for large super cells, 

the main determinant to get ripples is the applied strain. 
However, the structural patterns that occurs in defective 
graphene define a preferential direction and break the up- 
down symmetry that exits perpendicular to the layer. It 
is this preferential direction and the breaks in the sym- 
metry that are explained by the reconstructed vacancies, 
which break the hexagonal symmetry of graphene with 
their goggles-pentagon structure. Vacancies thus play a 
key role in determining the shape and symmetry of the 
global deformation patterns of graphene. 

We simultaneously characterize the local curvature at 
the vacancies and the height of atom 2 above the local 
tangent-plane in terms of the bending angle #12, as de- 
fined in [TJ In [2j>, we plot the bending angle for both 
defective and pristine graphene. Under moderate com- 
pression (1-2%), the local bending angle for defective 
graphene is about three times greater than that of pris- 
tine graphene and the graphene rippling is clearly made 
easier by the presence of the vacancies. The dependence 
of #12 on the size of the supercell, as shown in the inset 
of[2j>, further demonstrates the coupling between the lo- 
cal geometry of the defect and the global deformation, 
in that for 4 x 4 supercell, a strain larger than 1.0% is 
required to obtain an increase in #12 and to allow the 
corrugation of the layer to begin. However, for larger 
super cells, much smaller strains cause appreciable defor- 
mations around the vacancies. 

Energetics. The range of applied strain used here (a 
few percent) is comparable to that used in experiments 
in which an appreciable corrugation of graphene was re- 
ported for supported layers (36l Eol Elj , chemically func- 
tionalized graphene [16| , or defective layers [17j. As a 
result of the imposed periodicity and the finite-sizes used 
in our calculations, long- wavelength deformations appear 
primarily to be related to the strain, i.e. vacancies do not 
create significant corrugation in the relaxed geometry at 
the equilibrium lattice constant. However, we observe 
that for the range of strain considered here, the system is 
compressed and thereby assumes a rippled configuration 
at an energy between two and three times lower for de- 
fective graphene than for the pristine layer. For a 10x10 
supercell, the energy required to create ripples for a com- 
pression of 3% is reduced to almost half its value in the 
presence of vacancies, as seen in [3]A.. This difference in 
energy is consistent with the proposed role of vacancies 
as a source of ripples in recent experiments 

Magnetism versus strain: several spin solutions. In the 
foregoing part we commented on the fact that isotropic 
strain can be used to tune the vacancy structure and 
the curvature of the layer at the site of the defect. We 
now focus on the influence of these structural changes 
on the magnetic and electronic properties of the vacan- 
cies. We first consider the case of zero strain, i.e. at the 
equilibrium lattice constant. The usual reconstruction 
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FIG. 3. Variation in total energy per atom (A), spin moment 
(B), and bond length (C) versus strain for three different con- 
figurations of vacancy in the 10x10 supercell: high-spin flat 
(triangles), low-spin flat (squares) and low-spin rippled (cir- 
cles). The inset in panel (a) details the change in total energy 
with respect to the low-spin solutions at about zero strain. 
The filled and empty symbols in panel (C) represent the 1-2 
and 1-1 distances, respectively (see the inset of [T] for atomic 
labels). The marks given by stars refer to the average dis- 
tances as the geometry departs from the usual vacancy to 
the structure in[TjB. Note that the magnetism disappears at a 
strain of ^2% when allowing if out-of-plane deformations are 
allowed. 



[38|, [39|| is accompanied by a decrease in the 1-1 distance 
in [1] and an increase in the 1-2 distances. The two dan- 
gling bonds in the type 1 atoms are thus saturated while 
atom 2 remains uncoordinated. It is the polarization of 
the corresponding dangling bond that is the main reason 
behind the appearance of a spin moment of ^1.5/iB asso- 
ciated with the carbon monovacancy. However, we were 
able to stabilize another reconstruction characterized by 
a different structural distortion, in which the 1-2 distance 
decreases, while the 1-1 distance increases [see the local 
structures and distances in[3B,C]. This structure has a 
larger spin moment of ~1.82/i#. We refer to this latter 
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structure as the high-spin (HS) configuration, and to the 
former as the low-spin (LS) configuration. At zero strain 
the LS structure is more stable than the HS by around 
250 meV. 

The behavior of these two structures as a function of 
strain is summarized in [3j When applying tension to the 
layer, both the HS and LS structures remain flat and 
almost become degenerate. In both cases the spin mo- 
ments show a slight increase. Conversely, when the layer 
is compressed, the HS structure becomes unstable. In- 
deed, it is only possible to stabilize the HS configuration 
under compression provided that the layer is constrained 
to remain flat. Even for flat graphene, the energy differ- 
ence between the HS and LS states increases significantly 
with compression, and for strains greater than 1.5% the 
HS configuration spontaneously transforms into LS-flat 
(i.e. low-spin constrained to be flat). In both structures, 
the spin moments decrease as the layer is compressed. 
The change is greater for the LS-flat configuration, which 
varies from 1.98 \±b for a +3% deformation to 1.15 \±b at 
—3%. A more dramatic reduction in the spin moment is 
seen if ripples are allowed to form in the graphene layer. 
[3p shows that the spin moment of the LS-rippled (i.e. 
low-spin free to ripple) vacancy decreases sharply when 
the compressions exceeds 0.5%. In fact, the ground state 
of the monovacancy becomes non-magnetic for compres- 
sive strains in the range of 1.5-2%. 

The changes in the spin moments of vacancies may 
be understood by analyzing their electronic structure. 
When carbon atom 2 is restricted to remain in the plane, 
the spin moment is mainly associated with the sp 2 dan- 
gling bond, and the contribution from the p z states is 
smaller. The strain-induced deformation of the vacancy 
and subsequent out-of-plane displacement of atom 2 gives 
rise to the hybridization between the out-of-plane p z 
states and the in-plane sp 2 states. For flat graphene 
the spin-polarized impurity level associated with the va- 
cancy therefore has a strong sp 2 character and remains 
essentially decoupled from the delocalized electronic lev- 
els of graphene due to their different symmetries. For 
the rippled layer, these two types of electronic states are 
strongly hybridized. This hybridization results in the de- 
localization of the defect levels which eliminate the mag- 
netism and explains why the LS-rippled configuration has 



a lower spin moment. Expressed in three dimensions, the 
rippling is believed to transform the hybridization of the 
C vacancy atoms in sp 3 , thereby removing the cause of 
local magnetism, as seen at the higher compression of the 
case of[TJ>. 

Rippling of the graphene layer occurs in a range of 
cases. One such case is where graphene is deposited 
on substrates, in which there is a significant mismatch 
in between lattice parameters. For example, ripples 
have been previously observed for graphene deposited on 
Ru(0001) [36], as well as on other metallic substrates. 
Strain can also be applied and controlled by placing ex- 
foliated graphene on flexible substrates, as reported re- 
cently [42, 43]. Under such conditions, it may be possible 
to create defects in regions of different curvature using a 
focused electron beam, by applying a technique similar 
to that recently demonstrated by Rodriguez- Manzo and 
Banhart [44]. The existence and spatial distribution of 
the spin moment associated with the vacancies, and its 
variations with deformation, could then in principle be 
monitored using a spin-polarized scanning tunneling mi- 
croscope. [33, |45| 

In summary, we have shown that the magnetic and 
structural properties of carbon vacancies depend strongly 
on the local curvature induced by defects. Together with 
the global rippling pattern, this curvature can be con- 
trolled by the application of a strain. By compress- 
ing graphene by up to 2%, we can adjust the spin mo- 
ment of the vacancy to between and 1.5 \±b- At our 
high-compression limit, the rippling allows the defect C 
atoms to be arranged in a geometry that is different from 
the vacancy as generally observed. Such strains may be 
achieved experimentally [42|, |43[ and are comparable with 
those found when graphene is deposited on a range of 
different substrates. In fact, our results suggest that the 
magnetism that is induced in graphene by the presence of 
defects can be controlled using isotropic strain and other 
mechanical deformations. 
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